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The Charge Radius of the Proton 
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Recently, the charge radius of the proton was extracted for the first time from muonic hydrogen. 
The value obtained is five standard deviations away from similar measurements of regular hydrogen. 
This talk discusses work done in collaboration with Richard J. Hill, to address this discrepancy. First, 
we have studied the extraction of the charge radius of the proton from electron-proton scattering 
data in a model-independent way. We have shown that previous extractions, spanning a period of 
over 40 years, have underestimated their errors. Second, we have looked at a model-independent 
analysis of proton structure effects for hydrogen-like bound states, using the tool of an effective field 
theory, namely NRQED. We have identified hidden model-dependent assumptions in the theoretical 
I " calculation behind the muonic hydrogen result. 

o- 

; I. INTRODUCTION 

How big is the proton? To answer this question we first have to define what do we mean by the proton's size. 
The most common definition is of the charge radius of the proton. The matrix element of the electromagnetic 
\^ ' current between nucleon states gives rise to two form factors, 
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where q — pj—pi and m is the nucleon mass. F\ and F2 are often referred to as the Dirac and Pauli form factors, 
respectively. The Sachs electric and magnetic form factors are defined as Ge(q 2 ) — F\(q 2 ) + q 2 ^{q 2 ) / '4m 2 , 
Gi\i(q 2 ) = Fi(q 2 ) + F2(q 2 ). The values of Ge and Gm at q 2 = give the electric charge and magnetic moment 
'< of the nucleon, respectively. For example, G E (0) = 1 and G P M {Q) — \i p « 2.793. The electric charge radius of 
| [ the proton is defined by the slope of G P E 
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The motivation for this definition is that in the Breit, or brick-wall, frame, r p E is the root-mean-square radius 
of the charge distribution. Notice, though, that the definition ^ is independent of any specific frame. 

The most direct way to extract the charge radius of the proton is by measuring G P E from scattering exper- 
iments, usually electron-proton scattering. An indirect way is to use spectroscopy of a bound proton-lepton 
system. Intuitively, the fact that the proton is not a point particle but has a structure, modifies the potential 
seen by the leptons. For a point particle the amplitude for p+£ — > p+£ is proportional to 1/q 2 , leading to the well 
known —Za/r potential. Including q 2 corrections from proton structure lead to an amplitude proportional to a 
constant. In position space it gives a delta function type of a potential, proportional to the charge radius squared, 
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- . ■ AirZa.5 3 (r)(r E ) 2 /6. The corresponding correction to the energy level is AE r P = 2(Za) 4 m 3 (r p E ) 2 S eo /3n 3 , where 
mi is the lepton mass (I = e, m r = m^m p / (m£ + m p ), and m p is the proton mass. For regular and muonic hy- 
drogen m r Ri me.. Since m p 3> m e , proton structure effects for muonic hydrogen are greatly enhanced compared 
to regular hydrogen. 

For a long time it was anticipated that a measurement of the Lamb shift in muonic hydrogen would reduce 
the error by an order of magnitude compared to measurements from electron-proton scattering and regular 
hydrogen spectroscopy Last year, the CREMA collaboration at the Paul Scherrer Institute in Switzerland 
have reported the first measurement of the charge radius of the proton from the Lamb shift in muonic hydrogen 
0. The value obtained was r p E = 0.84184(67) fm. As expected, the error is an order of magnitude smaller 
than the value of r p E = 0.8768(69) fm, extracted mainly from (electronic) hydrogen spectroscopy by CODATA 
(Committee on Data for Science and Technology) ||. But the central value is five standard deviation away! 
It is quite rare in high energy physics to have a five sigma deviation, which motivated several studies of new 
physics effects [4j-l6|. How should we interpret this result? 
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II. MODEL INDEPENDENT EXTRACTION OF THE PROTON CHARGE RADIUS FROM 

ELECTRON SCATTERING 



Since we are having two conflicting measurements, it is useful to look at a third one, the measurement from 
electron-proton scattering data. What is the value extracted from that measurement? The first place to look is 
the particle data book. The 2010 edition [?J lists 12 extractions from electron-proton scattering data that span 
the period of 1963-2005 and the range of 0.8 — 0.9 fm, most with quoted uncertainties of 0.01-0.02 fm. These 
are based on different data sets and different extraction methods. It is hard to draw any conclusion from this 
list. Indeed, the Particle Data Group (PDG) itself does not use any of these extractions "...for averages, fits, 
limits, etc." 

Part of the problem is the use of different data sets. This problem can be overcome by a world effort in 
combining data, perhaps in a way similar to the Heavy Flavor Averaging Group Q . A more subtle issue is the 
question of how to extract the charge radius of the proton from the measured values of G v E (q 2 ). This seems 
like a simple problem. We have a data set, and we just need to find the slope at q 2 — 0. The problem is that 
we do not know the functional form of Gs(g 2 ), so the extrapolation to q 2 — is not as simple. Since we do 
not know the functional form, one would like to extract the charge radius in a model-independent way. This 
rules out some extractions based on a specific form of the form factor, typically a sum of poles and continuum 
states. What we would like to have is a general paramctrization, a series expansion for example, that gives a 
consistent result for the charge radius as we fit more parameters to the data. 

There are several series expansion that are used. Perhaps the simplest one is to expand GE{q 2 ) in a taylor 
series around q 2 = 0. Many of the sources mentioned by the PDG have used such an expansion, but typically 
with a small (2-3) number of parameters. Another type of a series expansion is the continued fraction expansion 
suggested in @. In 0, we have advocated the use the so called "^-expansion" , which is a standard tool in 
analyzing meson transition form factors, (for references see fiof). The z-expansion utilizes the known analytic 
properties of G p E (t), (t = q 2 ): it is analytic in the complex i-plane outside a cut of t e [4m 2 , oo]. The electron- 
proton scattering data lies on the real t < axis. It is therefore common to define Q 2 = —q 2 = —t. One can 
map the domain of analyticity onto the unit circle by the conformal transformation 
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where t cut — 4m 2 and to is a free parameter representing the point mapping onto z = 0, see figure [1] . It is 
convenient to choose to = 0, but physical observables are independent of the choice of to- The maximum size 
of z depends on the distance between i cut and t — 0. This implies that if we can increase i cu t, the maximum 
size of z decreases and the convergence of the series is improved. Since G P E is analytic in the unit disc of the z 
plane, it can be expanded as a taylor series in z: G E (q 2 ) = 53fc=o ak z (q 2 ) ■ 
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FIG. 1: Conformal mapping of the cut plane to the unit circle. 



Does it matter which expansion we use? This question was studied in [Toj . Using the low q 2 data tabulated 
11], we have compared the extraction of the charge radius of the proton using 4 different expansions: Taylor, 



continued fraction, z expansion, and z expansion with a bound on the coefficients, \ak\ < 10. While fits of the 
different expansion with 2 parameters agree well, as we increase the numbers of parameters, the errors for the 
first 3 fits grow without bound. Only the constrained z expansion is stable and therefore mo del- independent. 
Another interesting feature is that the error on the extracted charge radius for a fit with 2 parameters is actually 
an underestimate. The final answer has an error that can be larger by a factor of 2. The conclusion is that 
many previous extractions have underestimated their errors. 

Since the z expansion is related to the analytic structure of the form factor, information about the imaginary 
part of G E (t + i0) directly translates into information about the expansion coefficients. The analysis of [lCj had 
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shown that a bound of |ofc| < 10 is very conservative. Using the electron-proton scattering data tabulated in 
[l2|, we have found that r p E — 0.870 ± 0.023 ± 0.012 fm, where the first error is obtained using \ak\ < 5, and the 
additional error for a fit assuming |a/-| < 10. The errors might seem large compared to some of the values listed 
in Q, but as we have shown in [To], this can be attributed to the model-dependence of the previous extractions. 

The error can be reduced by including neutron and tttt data. The inclusion of neutron data allows for a 
separate fit of the isoscalar and isovector components of G P E . For the isoscalar component, the cut starts at 
t = 9m 2 . By including the irir data, namely irir production and ittt — > NN, we can effectively raise the isovector 
threshold to t = 16m 2 . As explained before, these typically lead to a better convergence and a smaller error. 
The best value obtained in 110] is r p w = 0.871 ± 0.009 ± 0.002 ± 0.002 fm, where the first and the second error 
are as above, and the final er7or is from the tttt continuum contribution, see for details. 

The problem of model-dependent extraction of such basic non-perturbative parameters is not unique to the 
vector form factor. The matrix element of the axial current between nucleon states gives rise to the axial form 
factor. Analogously to ([2]), its slope at q 2 = defines the nucleon axial radius r^, or alternatively, the axial 
mass defined as = y/l2r A . Extractions of the axial mass usually use a dipole model for the axial form 
factor, F A = F A (0) [1 - q 2 /(m A po ) 2 ]" 2 . Recently, MiniBooNE has extracted the axial mass from neutrino 
scattering data. The value obtained, m A ipole = 1.35 ± 0.17 GeV [l3j], is in conflict with the value obtained from 

pion electro-production data, m A polc = 1.07 ± 0.02 GeV [Hj]. Notice the superscript "dipole", emphasizing 
that fact that a dipole model was used. Using a model- independent approach [15|], we have obtained the value 
tua = 0.85^q'q 2 ± 0.09 GeV. From an illustrative dataset for pion electro-production, we have extracted the 
value, wia — 0.92^ \^ ± 0.08 GeV. As a check, we have also extracted m A polc from the same data sets by using 

the dipole model. We obtain to^' p ° 1c = 1.29 ± 0.05 GeV (neutrino scattering) and m A polc = 1.00 ± 0.02 GeV 
(pion electro-production). As far as the axial mass is concerned, the discrepancy can be attributed to the use 
of dipole model. 

Returning to the issue of the charge radius of the proton, the value obtained from a model-independent 
extraction from scattering data is r p E — 0.871 ± 0.009 ± 0.002 ± 0.002 fm. It is more consistent with the 
CODATA value, r p E = 0.8768(69) fm, than the CREMA value, r p E = 0.84184(67) fm. Let us look more carefully 
at the muonic hydrogen extraction. 



III. MODEL INDEPENDENT ANALYSIS OF PROTON STRUCTURE FOR HYDROGEN-LIKE 

BOUND STATES 

The CREMA collaboration has not measured the charge radius of the proton directly. They have measured 
the Lamb shift of AE = E[2P^ - E[2S£=%] = 206.2949 ± 0.0032 meV, between the 2£jfr, 1 2 ) and 2P]f 3 = / 2 2 ) 
energy levels. To extract the charge radius it should be compared to the theoretical expression [l(| H3] 

AE = 209.9779(49) - 5.2262(r|) 2 + 0.0347(r|) 3 meV , (4) 

to obtain r p E — 0.84184(67) fm. The constant term arises from many contributions, mostly QED effects related 
to the muon. We have explained in the introduction how the quadratic term arises. Where does the cubic term 
come from? 
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FIG. 2: Two photon exchange 

In the calculation that leads to (Q} , the cubic term arises from a two photon exchange between the lepton and 
the proton, see figure 2. To obtain the cubic term, it is separated to a "proton" and "non-proton" contributions. 
To calculate the proton part, one uses "Feynman rules" in which onshell form factors are inserted into the vertices 
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of Feynman diagrams for a relativistic pointlike particle. The proton part of the amplitude is now a function 
of Ge and Gm- Assuming a dipole model for the form factors d — Gj(0) [1 — q 2 /A 2 ]~ 2 , the amplitude can 
be expressed as function of one parameter, A. Since one parameter model for the form factors is assumed, this 
parameter can be related to r p E . "Expanding" in powers of r p E , an approximated cubic dependence is obtained. 
Not expanding in powers of r E , and using A 2 = 0.71 GeV, it leads to AE = 0.018 meV [l|. To explain the 
discrepancy it should be 0.258(90) meV, if using r p E from scattering data 10], or 0.311(63) meV, if using r v E 
from spectroscopy Q. 

This procedure raises several questions. First, can we separate the two photon amplitude into a "proton" 
and "non-proton" contributions? Second, why is the insertion of the form factors in vertices valid? Third, 
even if the insertion were valid, the dependence on the cubic power of the charge radius is only valid for a one 
parameter model of the form factors and in the large proton mass limit [l6j . Since we are discussing precision 
measurements, a more solid theoretical treatment is clearly needed. 

To study proton structure effects in hydrogen-like bound states in a model-independent way, one can use an 
effective field theory, namely Non Relativistic QED (NRQED) [18]. The NRQED Lagrangian up to order 1/m 3 

is fUEt 
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We also need the contact interaction 

ipl(Ttp p ■ i>Jatpi iplippipjipi 

i-coutact — «1 r 0,% ■ (O) 

miirip miirip 

The knowledge of the Wilson coefficients allows us to determine the proton structure corrections to the energy 
levels. Thus 

xjp, ^ x ml{Zaf ( Zom proton d 2 \ 

SE {n, i) = Sn-^p- {-^/n -—)■ (7) 

The Wilson coefficients are determined as follows. Coefficients of operators that contain coupling to one photon 
field, namely the Cj's in ([5]) apart from cai and CA2, are determined by the derivatives of the form factors, F\ 
and F2, at q 2 = 0. Coefficients of operators that contain coupling only to two photons, namely cai and CA2, are 
determined by the forward and backward Compton scattering. Coefficients of the contact terms, namely d\ and 
di-, are determined by using the zero- momentum limit for the p + £ — > p + £ amplitude. The tree level, 0(a), 
amplitude is reproduced by the effective field theory, and d\ and d 2 receive a non-zero contribution starting at 
0{a 2 ). Here we consider only the spin-independent coefficient, d 2 - The relevant proton matrix element is the 
forward Compton amplitude 

^ij d 4 xe^^k,s\T{J^ m ix)J^ m X0)}\k,s) = LgT + ^ W 1 + (V - (k» - ^jf-^J W 2 , 

(8) 

where the Wi's are functions of v = 2k ■ q and Q 2 = —q 2 . The matching condition is 
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where A is the photon mass and 5d 2 denotes the contribution to d 2 in addition to the point particle value. To 
determine dd 2 , we need to know W\ and W 2 . How can we determine them? 

The imaginary part of W\ and W 2 can be related to measured quantities. Inserting a complete set of states 
to ©, the imaginary part of the Wi's can be divided to proton and non-proton contribution. The proton 
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contribution can be expressed in terms of form factors, while the non-proton contribution can be expressed 
in terms of inelastic structure functions. To reconstruct the W^s from their imaginary parts, we need to use 
dispersion relations: 

w JvcMQ 2 ) 2 v' 2 (v' 2 -v-) 

1 f°° 
W 2 {v,Q 2 ) = - \ 



t (Q 2 



(10) 



For convergence, W\ requires a subtraction. As a result, one needs to know also Wi(0, Q 2 )- Currently IUi(0, Q 2 ) 
cannot be calculated from first principles or extracted from experiment. This fact introduces an uncertainty 
that was largely ignored in the literature. 

Can we say anything model-independently about Wi(0, Q 2 )7 We can calculate it in the large and small limit 
of Q 2 . For small Q 2 the photon "sees" the proton almost like an elementary particle. Using NRQED to calculate 
Wi(0, Q 2 ) up to 0(Q 2 ) (including) Q3 

O 2 

Wi(0, Q 2 ) = 2(4 - 1) + 2^-| (c Al +c 2 F - 2c F c wl + 2c M ) , (11) 

where we have described above how the Cj's are determined. For large Q 2 the photon "sees" the quarks inside 
the proton. Using the Operator Product Expansion to calculate Wi(0,Q 2 ), we find that Wi(0,Q 2 ) ~ 1/Q 2 
for large Q 2 [l8j |. The intermediate Q 2 region is not constrained by existing measurements. This introduces 
model-dependence into the theoretical prediction. 

We have described above the procedure leading to ((H). We can now address some of the questions that 
it raises. First, the separation to proton and non-proton states is only justified for the imaginary part of 
the Wi'a. Wi(0,Q 2 ) cannot be rigorously divided in such a way. Second, the "proton part" of Wi(0, Q 2 ) 
is evaluated in the literature by inserting onshell form factors into the vertices of Feynman diagrams for a 
relativistic pointlike particle. The result, which we refer to as the "Sticking In Form Factors" (SIFF) model, 
gives Wf IFF (0,Q 2 ) = 2F 2 (Q 2 ) [2F X (Q 2 ) + F 2 (Q 2 )]. This result is not derived from a well defined local field 
theory. In fact, no local Lagrangian can give such Feynman rules. Furthermore, the resulting expression for 
Wi(0, Q 2 ) has the wrong large Q 2 behavior. The "non-proton" part is obtained by multiplying part of the 
term proportional to cai in (|lip by a function of Q 2 (lq |. This function is not derived from first principles and 
again does not have the correct large Q 2 behavior. We stress again that unlike the imaginary part of the Wi's, 
the separation of W\(0,Q 2 ) into proton and non-proton parts is not well-defined. In conclusion, the standard 
calculation behind (HJ is model-dependent. One might ask though, how severe is the model dependence? 

Considering (flUl) . it is natural to decompose the two photon contribution to the energy levels as Ai?*™ -7 — 
A£; continuum + A £ P roton + ^gWi (o,Q 2 ) _ The firgt twQ termg ^ this equation can be extracted from data and are 

therefore model- independent. The last term is model-dependent. We look at the n = 2 level of muonic hydrogen. 
The continuum contribution was recently determined to be [22j], A£^™ tlnuum = 0.0127(5) meV. To illustrate 
the magnitude of the proton contribution, let us use a simple dipole model: Ge(q 2 ) ~ Gm(? 2 )/Gm(0) ~ 
[1 - q 2 /A 2 }- 2 with A 2 = 0.71 GeV 2 . We find AE p J° ton = -0.016 meV. This result can be refined by using the 

measured values of the form factors. In the SIFF model, using the dipole form factors, we find AE^^ '® '' SIFF = 
0.034 meV. The sum of the last two terms gives the 0.018 meV mentioned above |l|. In other words, the 
theoretical prediction for the "proton contribution" consists of a sum of two terms. The model dependent piece 
is twice as large as the model-independent piece and has the opposite sign. Needless to say, this result is very 
alarming. The theoretical prediction used for the extraction of the charge radius of the proton from muonic 
hydrogen relays heavily on the validity of the SIFF model. It is not hard to construct models of Wi(Q, Q 2 ) that 
would that have the correct small-Q 2 and large-Q 2 behavior, but give a much larger contribution than the SIFF 
model. 



IV. CONCLUSIONS 



In this talk we have discussed the recent discrepancy between the extraction the proton charge radius from 
muonic {r p E — 0.84184(67) fm) and regular {r p E — 0.8768(69) fm) hydrogen. We have looked at the extraction of 
the charge radius from electron-proton scattering data. Such an extraction is complicated by the fact that the 
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functional form of the form factor is unknown. We have established in .10] that previous extractions suffer from 
model-dependence and that the error on the charge radius was underestimated. This situation is not unique to 
the vector form factor. The recent discrepancy in the extraction of the axial mass from the MiniBooNE data 
is in conflict with previous pion electro-production data, when a one-parameter model for the axial form-factor 
is used. When a model independent approach is taken, [l5j . there is no discrepancy, as far as the axial mass 
is concerned. We have used a model independent approach to extract the charge radius of the proton from 
proton, neutron and irir data. The value of the charge radius we obtain, r p E = 0.871(11) fm, is closer to the 
regular hydrogen result. This motivates a closer look at the muonic hydrogen determination. 

We have presented the NRQED formalism for isolating and controlling proton structure effects in hydrogen- 
like bound states. The most severe problem with the theory behind the muonic hydrogen measurement arises 
from the contact interaction © . The Wilson coefficient d 2 is not determined by proton form factors or inelastic 
structure functions. To determine it, the two photon amplitude is needed. 

Let us summarizes the discussion of the two photon amplitude. To determine di we need the imaginary part 
of the W^s which can be extracted from data, and also Wi(0,Q 2 ) which currently cannot be extracted from 
data. Unlike the imaginary part, Wi(0,Q 2 ) cannot be written model independently as a sum of "proton" and 
"non-proton" terms. The theoretical calculation behind the muonic hydrogen result assumes such a separation 
and relies on the SIFF model to estimate the "proton" part. Even worse, the SIFF model is not identified 
as a model, but presented as a theory. In [l8| we have calculated model-independent properties of W\(Q,Q 2 ), 
namely its low Q 2 behavior via NRQED and high Q 2 behavior via OPE. The intermediate region is poorly 
constrained. The lack of theoretical control over W\{Q,Q 2 ) introduces theoretical uncertainties that were not 
taken into account in the literature. 

While not resolving the discrepancy per se, our studies have exposed hidden model-dependence in the ex- 
traction of the charge radius of the proton from scattering data and in the theoretical calculation behind the 
extraction of the charge radius of the proton from muonic hydrogen. Hopefully, these studies will serve as a 
basis for a more careful extraction of this basic quantity. 
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